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ABSTRACT 

Discs that contain dead zones are subject to the Gravo-Magneto (GM) instabihty that 
arises when the turbulence shifts from gravitational to magnetic. We have previously 
described this instability through a local analysis at some radius in the disc in terms 
of a limit cycle. A disc may be locally unstable over a radial interval. In this paper, 
we consider how the local instability model can describe global disc outbursts. The 
outburst is triggered near the middle of the range of locally unstable radii. The sudden 
increase in turbulence within high surface density material causes a snow plough of 
density that propagates both inwards and outwards. All radii inside of the trigger 
radius become unstable, as well as locally unstable radii outside the trigger radius. In 
addition, a locally stable region outside of the trigger radius may also become unstable 
as the gravitational instability is enhanced by the snow plough. For the circumstellar 
disc model we consider, we find that a quarter of the disc mass is accreted on to 
the central object during the outburst. The radius out to which the disc is globally 
unstable is twice that for which it is locally unstable. 

Key words: accretion, accretion disks - planets and satellites: formation - proto- 
planetary disks - stars: pre-main sequence 



1 INTRODUCTION 

Accretion discs are ubiquitous in the Universe, they form on 
all scales from planetary, to stellar, to galactic. They trans- 
port angular momentum outwards allowing material to spi- 
ral inwards through the disc (e.g. |Pringlelll98ll ). If the disc 
is fully ionised, the magneto rotational instability (MRI) 
drives turbulence and t hus angular momentum transport 
(JBalbus fc Hawlevlll99J ). However, below a critical level of 
ionisation, a dead zone may form at the mid-plane where the 
MRI does not operate |Gammiel 119961 : iGammie &: Menoul 
1 19981 : [Turner fc Sandl2008l '). The inner parts of the disc are 
hot enough to be thermally ionised. However, further out, 
only the surface layers of the disc may be ionised by exter- 
nal s ources such as cosmic ray s or X -rays from the central 
star (JGlassgold. Naiita fc Igeal |2004 ). The restricted flow 
through the disc causes a build up of material in the dead 
zone. If there is sufficient accretion infall, some turbulence in 
the dead zone may be driv en by self-gravity when the layer 
becom es massive enough (iPaczvnskil 119781 : iLodato fc Ricd 
[2OOJ). However, if it does not transport angular momentum 
fast enough, the gravo-magneto (GM) disc instability can 
result. This occurs when the temperature in the dead zone 
reaches the critical required for t he MRI to be triggered and 
an accretion outburst ensues (|Armitage. Livio fc Pringld 
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2OOII: IZhu. Hartmann fc Gammid |2009| . l2010al : IZhu et all 
2010bl : iMartin fc Lubowll2011al '). 



The GM disc instability can be explained as transitions 
between steady state disc solutions plotted on a state dia- 
gram showing the accretion rate t hrough the disc against th e 
surface density for a fixed radius (jMartin fc Lubowll2011al V 
There are three steady state disc solutions, two fully MRI 
turbulent solutions (one that is thermally ionised and one 
that is fully ionised by external sources such as cosmic rays 
or X-rays) and the GM solution. The GM solution consists of 
a self gravitating, MRI inactive mid-plane region with MRI 
active surface layers such that the flow through the disc is 
at the steady rate. For a given radius, there may be a range 
of accretion rates for which no steady solution exists. If the 
accretion rate on to the disc lies in this range then the disc 
will be locally GM unstable. The disc transitions between 
the steady thermal MRI and GM solutions cause repeat- 
ing outbursts over time. This is sim ilar to the "S-curve" 
used to explain dwarf nova outbur sts (jBath fc Pringlg|l982l : 
iFaulkner. Lin fc Papaloizoulll983l 'l. 

For a range of radii in the disc, there exists a local limit 
cycle. I n this work, we extend the model of lMartin fc Lubowl 
(|2011al 'l to consider how the local limit cycles at different 
radii are coordinated and thus how the global disc evolu- 
tion operates. We examine the radius at which the outburst 
is triggered and how far the outburst propagates through 
the disc. We flnd that the radii are connected through 
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a snow plough effect where the dense material is pushed 
both inwards and outwards during the outburst. Previous 
works have considered the propagation of the outb urst (e.g. 
lArmitage. Livio fc Pringi3l200ll : IZhu et al.ll2010bl ). but not 
within the framework of the limit cycle. 

The GM disc instability may occur in accretion 
discs on many scales. It is thought to describe FU 
Orion is outbursts in discs around yo u ng stellar objects 
(e.g. lArm itage. Livio fc Pringld [JOOll: IZhu et al.l l2010bl : 
iMartin et al. 2012allbl: Istramatellos et al.ll2012h . It is also 
likely to occur in circumplanetary discs that form around 
massive planets as they form within the circums tellar disc 
l|Martin fc Lubowll2011]bl : lLubow fc Martinll2012al lbl). In this 
work we restrict our analysis to that of a circumstellar disc, 
but note that the underlying mechanism is the same for all 
scales. 



2 LOCAL DISC INSTABILITY 

The layered disc mod el we apply is based on 
lArmitage. Livio fc Pringld (|200ll) and further deve l- 
oped in IZhu et al.l (I2OI0S) ar id [Martin fc Lubowl (|2011al ). 
Following iMartin fc Lubowl (|2011al ). we first solve the 
steady state layered accretion disc equation 



M = 37r(l/mEm + I'gSg) 

with the steady energy equation 
9 

aT = -^ (rml^niSm +TI^gEg 

o 



(1) 



(2) 



where M is the steady accretion rate through the disc, T 
is the midplane temperature, Eg is the surface density in 
the dead zone layer, Em is the MRI active surface density, 
the total surface density is E = Em + Eg, fi is the Keplerian 
angular velocity, Tm is the optical depth of the magnetic layer 
and T is the optical depth of the whole disc. The viscosity 
due to the MRI is 



Q. 



(3) 



where Cm is t he sound speed in the acti ve surface layers 
and Qm is the IShakura fc SunvaevI (|l973l ) a parameter in 
the layer. The disc is self-gravitating if the iToomre (1964) 
parameter is less than the critical, Q < Qait ~ 2, and in 
this case fg depends on Q. The viscosity due to self-gravity 
is 



^g - "g j^ , 
where Cg is the sound speed in the mid-plane layer and 



(4) 



Qcrit 



(5) 



for Q < Qcrit and zero otherwise JLin fc Pringldll987l . ll990h . 
We find that the exact form of a^ does not affect the disc 
evol ution provided it is a steeply declining function of Q (see 
also lZhu et al.ll2010bl . and the Discussion section). 

In this work we consider a constant infall accretion rate, 
Afinfaii, on to the disc. For a steady disc solution, the infall 
accretion rate is equal to the accretion rate through the disc, 
M, at all radii. However, if at some radius there is no steady 
solution, then the disc is unstable to the GM disc instability. 



As described in the Introduction, there are three differ- 
ent steady solutions, two MRI active solutions and the GM 
solution. The fully MRI active solutions are found by set- 
ting I'g = 0. They exist in two regions of the disc, first in the 
inner parts where the mid-plane temperature is sufficiently 
high for thermal ionisation (T > Tcrit, where Tcrit is the 
temperature above which the MRI operates) and secondly 
in the outer parts where the surface density is sufficiently 
small that cosmic rays or X-rays penetrate the whole disc 
(E < Ecrit, where Ecrit is the critical surface density that is 
ionised through external sources). The third steady solution 
is the GM solution, where there is gravitational turbulence 
near the mid-plane and MRI turbulence in disc surface lay- 
ers. This solution occurs where the mid-plane temperature 
is too low for the MRI (T < Tcrit) and the surface density 
is greater than the critical that can be ionised by external 
sources (E > Ecrit), and large enough for the disc to be 
self-gravitating. The steady state solutions at a given radius 
cover a range of accretion rates. In general, there are cases 
where there is a gap in accretion rates for steady state so- 
lutions. Accretion rat es that fall in that gap a t some radius 
are locally unstable JMartin fc Lubowll2011al ). For a fixed 
accretion rate on to the disc, there is generally a range of 
radii where that rate falls within the local gaps and thus a 
range of radii where the disc can be locally unstable. 

In this work we take an example of a circumstellar 
disc, but note that the same mechanisms apply to circum- 
planetary discs on a smaller scale. We choose a star of 
mass M = 1 M© and a disc with a IShakura fc SunvaevI 
(1973 ) viscosity parameter Om = 0.01, Tcrit = 800 K and 
Ecrit = 200gcm~^. Fig. [1] shows the state diagrams of the 
steady solutions at various radii in the disc. The surface tem- 
perature is shown against the surface density for the steady 
state solutions at a fixed radius. The surface temperature in 
the disc is related to an effective, steady accretion rate, Ms, 
with 



Te 



8-Kcr 



(6) 



All radii shown, except R = 0.1 au, show the three 
steady branches. There are two fully MRI active solutions, 
the upper branch with T > Tcrit and the lower branch with 
E < Ecrit . The third branch is the self gravitating GM solu- 
tion. The steady state solutions are shown in the solid thick 
lines and the dead zone branch (that is not a steady branch) 
is shown as the dashed line. The branches are labeled in the 
state diagram at _R = 5 au. The shaded region is where there 
is no steady state solution and the disc is locally unstable. If 
the imposed disc accretion rate (infall accretion rate onto the 
disc), represented as the dotted line, lies within the shaded 
region, then the disc is locally unstable at that radius. The 
state diagram at a radius of _R = 0.1 au does not have an un- 
stable region because a fully MRI active disc solution exists 
at all accretion rates because either E < Ecrit or T > Tcrit. 



2.1 Innermost Unstable Radius 

For a disc with an unstable region, the innermost radius that 
is locally unstable occurs where the thermally ionised fully 
MRI active solution has a mid-plane temperature equal to 
the critical temperature. The surface density is larger than 



© 2013 RAS, MNRAS 000, [lH8] 



Propagation of the Gravo-Magneto Disc Instability 3 



100 



100 




1000 lO'' 

2/g cm-2 



100 




100 



1000 10" 

2/g cm"^ 



10= 



1000 



100 r 



100 



Dead Zone 
Lower MR! 

I I I 



100 



1000 lO" 

2/g cm-2 



1000 10* 

2/g cm"^ 




Figure 1. The state diagrams at radii R = 0.1, 1, 2, 3, 5 and 10 au. The thick soUd Unes show the analytical steady states. The thick 
dashed line shows the dead zone branch. The thin lines with arrows show the global time-dependent numerical models. The dotted line 
shows the infall accretion rate that is equal to the steady state accretion rate through the disc. The shaded region is where there is no 
steady state solution and the disc is locally unstable. The star shows where the outburst is triggered globally and the cross where it is 
triggered locally. 



© 2013 RAS, MNRAS 000, [iHl] 



R. G. Martin & S. H. Lubow 



the critical that is ionised by external sources, if the unstable 
region exists. 

This location can been seen in the state diagrams to 
occur at a radius where the accretion rate through the disc 
is such that it lies on the lower end of the upper MRI branch 
(similar to the middle left plot in Fig.[T}. That is, the radius 
where the dotted line intersects the lower end of the upper, 
solid sloped line. The upper MRI branch begins where the 
mid-plane temperature is equal to the critical, T = Tcrit. 
We scale the variables to M' = M/(10"'' Mq yr"^), M' = 
M/Mq, a' = a/0.01 and T^^ = Tcrit /800 K. We find steady 
state fully turbulent solutions (with Em — E and v^ — 0). 
The radius at which the fully turbulent solution has the 
critical mid-plane temperature is 



R = 1.87 Trnrr au. 



(7) 



This is the innermost radius that is locally unstable to the 
GM instability. 

2.2 Outermost Unstable Radius 

The outermost radius for which no steady solution exists 
(and thus the disc is locally unstable to the GM instability) 
occurs where the GM solution has a mid-plane tempera- 
ture equal to the critical. In the state diagram this is where 
the accretion rate through the disc lies at the upper end 
of the GM branch (similar to the lower left plot in Fig. [T]). 
This is where the self gravitating, GM solution reaches the 
critical mid-plane temperature. We again solve the steady 
state equations, but now include the self gravitating term. 
However, the equations must be solved numerically. For an 
accretion rate of Minfaii = 10~^ M© yr~^, we find the radius 
to be 7? = 4.51 au. Thus, the locally unstable region for this 
disc extends from R = 1.87 au to 7? = 4.51 au. Outside of 
this region there is a steady GM disc solution, and inside 
there is a fully MRI solution. We note that the innermost 
and outmost unstable radii in the disc do not necessarily 
correspond to the inner and outer edges of a dead zone in a 
time-dependent disc. 

2.3 Outer Transition to MRI Branch 

The outer parts of the disc are fully MRI active where the 
surface density is smaller than the critical, E < Ecrit. This 
occurs at a radius 

R - 29.15 ^,^^^^ au. (8) 

Hence, in summary, in the region R < 1.87 au the disc has 
a fully MRI active solution, 1.87 < _R/au < 4.51 the disc 
has no steady solution (and is locally unstable to the GM 
instability). In 4.51 < _R/au < 29.2 there is a steady GM 
solution. In _R > 29.2 au the disc is fully MRI active. In 
Figs. [2] and 3 we label these regions as a function of radius 
in the disc and shade the region for which there is no steady 
solution. 



3 GLOBAL DISC INSTABILITY 

We solve the time-dependent accretion disc e quations tha t 
consist of a surface density evolution equation (|Pringlelll98lh 



and an energy equation (|Pringlg| 19861 : ICannizzd 1993lV The 
equatio ns are described in more detail in lMartin fc Lubowl 
(|2011al 'l (see their equations 1 and 3). We take a grid of 
200 points equally spaced in log R that extends from Rjn = 
5Rb to -Rout = 40 au (e.g. lArmitage. Livio fc Pringlj|200ll : 
iMartin et alj|2007l) . At the inner boundary we have a zero 
torque boundary condition to allow the inward fiow of gas on 
to the star. At the outer boundary, we choose a zero radial 
velocity condition to prevent material from leaving the disc 
there. There is a constant infall accretion rate on to the disc 
of Afinfaii ~ 10~^ M0 yr~^ at a radius of 35 au. As shown in 
the previous section, for this accretion rate the disc is locally 
unstable for a range of radii. 

The disc evolves through several outbursts before the 
limit cycle repeats. Thus, the initial conditions of the disc 
have no effect on the limit cycles shown. Fig. [l] shows the 
local state diagrams, for various radii, of the surface tem- 
perature against the surface density. The thick lines show 
the analytic steady state solutions (as described in the pre- 
vious section) and the thin lines show the numerical time- 
dependent evolution. The dashed line shows the steady infall 
accretion rate. 

The crosses mark the position at the time when the 
outburst is triggered at that radius (i.e. where T > Tcrit), 
and the stars mark the position when the outburst is first 
triggered globally. For the i? = 3 au case, the two positions 
coincide with the end of the GM branch and hence also with 
a temperature of T = Tcrit- However, at other radii, the posi- 
tion when the outburst is triggered locally does not coincide 
with the position where it is triggered globally. Once the 
outburst is triggered somewhere in the disc, it propagates 
very quickly and adjacent radii undergo their limit cycles 
similarly fast. At _R = 0.1 au, the mid-plane temperature of 
the disc does not drop below the critical during the cycle, 
hence there is no cross shown. Similarly, at i? = 10 au, the 
mid-plane temperature during the cycle does not exceed the 
critical temperature and so there is no cross shown. 

The outburst propagates inwards and hence all radii 
inside of the trigger radius show outburst behaviour. The 
total mass of the disc before the outburst is 0.45 Mq of 
which 0.34 Mq is within the dead zone and 0.024 Mq is not 
self-gravitating. During the outburst, a mass of 0.1 Mq is ac- 
creted on to the star. This amount of material is initially dis- 
tributed up to a radius of 2.9 au (that roughly coincides with 
the trigger radius). The outburst also propagates outwards. 
The dead zone has a high surface density that suddenly be- 
comes more turbulent when the MRI is triggered, causing 
material to move outwards. For example, at _R = 5 au, the 
disc has a steady state GM solution. However, this solution 
is not found because of the outward propagation. The disc is 
depleted after an outburst, leaving a less massive and cooler 
disc than the steady state would predict. The disc moves 
towards the steady state but the outburst is first triggered 
at an accretion rate lower than the steady rate. Even fur- 
ther out, where the steady solution is lower down the GM 
branch, the outward propagation is not sufficient to trigger 
an outburst and the disc remains close to the GM branch 
and then moves downwards after the outburst. The outer- 
most radius for which the disc moves up to the upper MRI 
branch during the outburst is around R — 9au, twice the 
radius of outer edge of the locally GM unstable region. 

Fig. [2] shows the disc structure leading up to an out- 
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Figure 2. Disc structure leading up to an outburst. Left: The total surface density. Right: The mid-plane temperature. The shaded 
regions show the unstable regions where no steady state disc solution exists. The regions are labelled where steady state, MRI and GM 
solutions exist. The profiles, in order of increasing height within the unstable region, are at times t = 1.22 X 10^ yr, 5.21 X 10* yr and 
1.12 X 10^ yr since the previous outburst. The final time shows the disc just before the outburst is triggered. The dashed line in the 
temperature plot shows the critical temperature, T^rit, above which the MRI drives turbulence. 



burst. The branch labels in the figure (Upper MRI etc.) 
refer to where steady state solutions exist. For the time- 
dependent evolution being plotted, these labels may not re- 
flect the instantaneous structure at each radius. The shaded 
region shows where there is no steady state disc solution. 
The innermost parts of the disc are fully MRI turbulent. 
However, beyond radius of R = 0.41 au, the mid-plane tem- 
perature of the disc drops below the critical value. At this 
radius, R — 0.41 au, there is a sharp increase in the surface 
density and this marks the inner edge of the dead zone. The 
narrow density peak is due to a local pileup of gas. It is 
the result of the outward spread of turbulent gas from the 
active, inner (inside 0.41 AU) region into the nonturbulent 
dead zone where it cannot spread outward any further. The 
narrowness is a consequence of the sharp transition criteria 
we apply between active and dead zones. If we applied a 
smoother transition criteria, the peak would be spread out 
over a larger range of radii. The dead zone covers some of 
the region where there exists a steady state MRI disc so- 
lution (as labeled Upper MRI). This overlap exists because 
the outer parts of the flow near _R = 1 au are limited to the 
accretion rate of the active surface layer. There is not suffi- 
cient flow in this region to reach the steady solution at the 
infall accretion rate. 

The local peak in the surface density within the unsta- 
ble region marks the radius in the disc at which the dead 
zone becomes self gravitating. This local peak moves inwards 
in time as material builds up. There is a similar local peak in 
the temperature profile that moves inwards slightly behind 
the surface density peak in time. Once this peak reaches the 
critical temperature required for the MRI, the outburst is 
triggered (as in the last time shown in Fig. [2|. The outburst 
is first triggered at a radius of 7? = 2.6 au. 

It is the outer parts of the dead zone that become self 
gravitating first (as shown in Fig. [2]). Hence, in the limit 



cycles in Fig. [T] the larger the radius, the sooner the disc 
reaches the GM branch. However, the speed at which the 
disc moves up this branch, decreases with radius. As a func- 
tion of radius, there is then opposite behavior in the time 
that the disc reaches this branch and the speed at which 
it travels up the branch. Hence, the outburst is triggered 
somewhere in the middle of the locally unstable region. 

Fig. [3] shows the radius dependent accretion rate 
through the disc at various times during the outburst cy- 
cle. The upper portion of each plot represents the inward 
flow and the lower portion the outward flow. The top left 
plot shows the accretion rate just before an outburst, that 
corresponds to the upper lines for the surface density and 
temperature shown in Fig. [21 There are three main parts to 
the disc, the inner part that is fully MRI active, the middle 
part where there is a dead zone and the large outer part that 
is self gravitating. In the very outermost parts of the disc, 
the accretion drops off because there is a zero flow outer 
boundary condition. This only affects the parts of the disc 
outside of the radius where material is injected, at 35 au. 

The top right hand plot shows the accretion rate shortly 
after the outburst is triggered. There are both inward and 
outward propagating waves moving away from the trigger 
point at 2.6 au. This is the snow plough effect, material 
moves both inwards and outwards away from the trigger 
radius. It is this snow plough that causes the various radii 
to be coordinated. The middle left plot shows the inner wave 
has propagated all the way in and the disc is now in outburst 
phase with a high accretion rate on to the star. The outward 
propagating waves continue to move outwards. In the mid- 
dle right hand plot, the dead zone has begun to re-form. 
This initially occurs at a radius of 8.1 au. The dead zone 
formation propagates both inwards and outwards as seen 
in the lower left plot. As the dead zone forms, the viscous 
torque there drops. Thus, material just interior to the dead 
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Figure 3. The radius dependent accretion rate through the disc during the outburst cycle. The plots shown arc at times t = 0, 140, 400, 
800, 1000 and 1400 yr from the start of the outburst. In each plot, the upper portion shows the inward accretion and the lower portion 
shows the outward accretion. 
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zone initially moves outwards because the outward viscous 
torque is larger than the inward torque. This is the cause of 
the narrow peak in the surface density at the inner edge of 
the dead zone (as described in Fig. 2). Eventually the inner 
parts of the disc move back to the quiescent phase as the 
dead zone blocks the flow, as seen in the bottom right hand 
plot. This marks the end of the outburst, although there is 
still some ongoing outward propagation in the outer regions. 
The material in the disc builds up again as shown in Fig. [5] 
until the cycle repeats. 



final GM disc outburst. If planets form between outbursts, 
they likely accrete on to the star during the outburst. As the 
infall accretion rate drops over time, the outbursts cease but 
a dead zone may remain. This dead zone provides a quies- 
cent region suitable for planet formation. The presence of a 
dead zone also affects the migration of planets through the 
disc. Planets embedded in low viscosity disc s can undergo 
very slow and sometimes chaotic migration (|Li et al.l l2009l : 
IYu et al.ll201Cl ). The nature of gap opening and migration 
of planets through a layered disc should be investigated in 
future work. 



4 DISCUSSION 

There are a number of simplifications in our model that 
should be investigated in fu ture work. Some of these were 
discussed in mor e detail in iMartin fc Lubowl (|2011al ) and 
iLubow fc Martini (|2012al V For example, we have assumed 
that the critical surface density that is ionised by external 
sources in constant with radius. However, a more realistic 
way to find the extent of the dead zone may be with a critical 
magnetic Reynolds number (e.g. Fle ming. Stone fc HawlevI 
I2OOOI : iMatsumura fc Pudritz 2003). The criti c al valu e how- 
ever remains uncertain (e.g. iMartin et al.l l2012a) . The 
ionisation may be further suppressed by effects such as 
ambipolar diffusion and the prese nce of dust and poly- 
cyclic aromatic hydrocarbons (e.g. Bai fc Goodman 200ff; 
Bai 2011; Pcrcz-Bcckcr fc Chiang 2011; Simon ct al. 2012; 
Dzvurkevich et al.l 120131 ). No matter how the dead zone is 



determined, the outburst mechanism works in the same way. 
The shape of the li mit cycle may be slightly different (see 
IMartin et alj|2012ah . but the principles remain the same. 

There remains some uncertainty in the prescriptions for 
the viscosities in the disc. Observations of the fully MRI 
turbulent discs i n dwarf novae and X-ray bi naries suggest 
Q f» 0,1-0.4 (|King. Pringle fc Livid 1200 it ). However, in 
discs around T Tauri stars, such a high a would cause th e 
outer discs to expand too quickly (Hartmann et al.lll99a ). 
Moreover, theoretica l work suggests a smaller a ~ 0.01 (e.g. 
iFromang et al.ll2007l ). We have taken Om = 0.01 in this work, 
but note that the qualitative evolution of the disc would be 
similar had we taken a larger value for dm- There is also some 
uncertainty in the value of Qg in a self gravitating disc. But 
its exact value has litt le influence of the disc struc t ure an d 
evolution, as noted bv lZhu. Hartmann fc Gammid (|2010al ). 
Within the self gravitating region the Toomre parameter in 
our model is nearly constant, Q « Qcrit- In our prescription 
the value of Og is not flxed and depends on the Toomre 
parameter, Q. In steady state, implicitly it varies with the 
local cooling timescale, icooi, such that ag oc l/(fcooif2) (see 
also Gammic 2001). 

In a real disc the infall accretion rate is not likely to 
be constant in time, but rather decrease. The locally unsta- 
ble region would decrease in time until the accretion rate 
is so low that a fully MRI active solution is found. For the 
parameters chosen in this work, this would be for accretion 
rates Mnfaii < 10"* Mq yr'^ 

The existence of dead zones in protoplanetary 
discs may have important implications for planet for- 
mation. Planets that survive af ter the disc has dis- 
persed (e.g. by photoevaportation iHoUenbach et al.l 1 1994 : 
lAlexander. Clarke fc Pringld I2OO6I ) likely formed after the 



5 CONCLUSIONS 

The GM disc instability can be explained as transitions be- 
tween steady state solutions in a state diagram that plots 
the surface temperature against the surface density at a 
fixed radius. We have examined local cycles at different 
radii. There is a range of radii for which a given disc is lo- 
cally GM unstable. For the parameters chosen in this work 



(Ecrit = 200 gcm-^ Tc, 



800 K, Minfaii = 10"" Mq yr" 



and a = 0.01), we find the unstable range to be 1.87 < 
R/a.u < 4.51. 

The outburst is first triggered in the middle of the range 
of locally unstable radii. For the parameters we chose, the 
trigger radius is at 2.6 au. The resulting snow plough of ma- 
terial that propagates both inwards and outwards links the 
different radii in the disc. All radii inside of the trigger radius 
become unstable and the mass in this region moves inwards 
to be accreted during the outburst. The snow plough effect 
pushes material outwards through the disc to radii which 
would otherwise lie along the steady GM branch and be lo- 
cally stable. The higher density causes these outer radii to 
also become unstable. The outburst propagates to a radius 
of around 9 au, about double that of the outermost locally 
unstable radius. Even further out than this, there are in- 
creases in density but these are not sufficient for a change 
in the state. 
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